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We report on the engineering of a non-dispersive (flat) energy band in a geometrically frustrated 
lattice of micro-pillar optical cavities. By taking advantage of the non-hermitian nature of our 
system, we achieve bosonic condensation of exciton-polaritons into the flat band. Due to the infinite 
effective mass in such band, the condensate is highly sensitive to disorder and fragments into localized 
modes reflecting the elementary eigenstates produced by geometric frustration. This realization 
offers a novel approach to studying coherent phases of light and matter under the controlled interplay 
of frustration, interactions and dissipation. 


Flat energy bands occur in a variety of condensed mat¬ 
ter systems, from the Landau levels of an electron gas, 
edge states of graphene pQ, Aharonov-Bohm cages in 
metal networks [2], frustrated magnets [3] to unconven¬ 
tional superconductors [4]. The common feature of these 
materials is the appearance of a divergence in the den¬ 
sity of states at the energy of the flat band, which pre¬ 
vents straightforward ordering. Consequently, any small 
perturbation may have a dramatic effect on the system. 
For example, interactions often lead to strongly corre¬ 
lated and exotic phases of matter as observed in the frac¬ 
tional quantum Hall effect 0, spin liquids [6 or spin 
ices Another important class of emerging phe¬ 

nomena in flat band systems originates from the effects 
of disorder, which are enhanced by the very large mass 
[9] and can significantly deviate from conventional An¬ 
derson localization. Examples include the inverse Ander¬ 
son transition (delocalization transition) [ID] , localization 
with unconventional critical exponents and multi-fractal 
behavior nn> and mobility edges with algebraic singular¬ 
ities [T2] . 

The observation of these phenomena in solid-state sys¬ 
tems is often complicated by extrinsic material-specific 
perturbations and the impracticality of engineering suit¬ 
able lattice geometries. Artificial lattices, recently imple¬ 
mented in a number of physical systems, allow simulat¬ 
ing this flat band physics in a controllable manner. Pio¬ 
neering works in photonic systems HM9] and cold atom 
gases |20if22] highlighted the key role played by geomet¬ 
ric frustration m and evidenced characteristic features 
such as the absence of wavepacket diffraction in a flat 
band [T7HT9] . But despite many predictions [231126] the 
specific role played by interactions or disorder has not 
been experimentally addressed in these artificial lattices 
so far. Also, steady-state bosonic condensation has not 
been realized in a flat band, and the question of the exper¬ 
imental properties of such condensate remains open. In 
the case of ultracold atomic systems such study remains 


delicate as flat bands usually appear at energies far above 
the ground state [ini na ng. Investigation of flat band 
condensation thus require challenging experimental tech¬ 
niques to implement complex-valued tunneling constants 
[2T1 or the coherent transfer of the atomic condensate into 
an excited state [22] . 



FIG. 1. a, Scanning electron micrograph of a ID Lieb (or 
Stub) lattice of coupled micropillars etched out of a cavity. 
The close up shows schematically the pillars structure, with 
two distributed Bragg reflectors (DBR) forming the cavity, 
and quantum wells (QW) inserted at the antinodes of the op¬ 
tical field, b, Geometric structure of the lattice: it contains 
three types of pillars (A,B,C) linked by couplings t\ , £2 to near¬ 
est neighbors and t' to next-nearest neighbors. An example 
of plaquette-state (see text) is highlighted in grey, c, Band 
structure calculated in the tight-binding model with t' — 0, at 
zero and finite spectral detuning between A and C sites (solid 
and dotted lines, respectively), d, Experimental real space 
emission at the flat band energy (conditions of Fig. 2Jd). e, 
Calculated density profile of a typical flat band mode (only 
one unit cell is shown). 


In the present work, we use exciton-polaritons to in¬ 
vestigate bosonic condensation in a flat band. These 
quasiparticles arise from the strong coupling between ex- 
citons confined in quantum wells and photons confined 
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in a semiconductor microcavity [29] . Their mixed light- 
matter nature allows efficient band structure engineering 
through their photonic component [30], while providing 
scattering channels and non-linearities through their exci- 
tonic component m- These assets have allowed studying 
polaritons in staggered [32, 33] . squared [34], honeycomb 
[16] and Kagome lattices [13] . In addition, their dissipa¬ 
tive nature makes them an archetypal non-hermitian sys¬ 
tem [35lf37] featuring novel dynamical universality classes 
[38] . In this Letter, we engineer geometric frustration by 
patterning a cavity into a ID Lieb lattice (also known 
as Stub lattice) [39 of optical micro-pillars. Photolumi¬ 
nescence at low excitation power directly reveals the pres¬ 
ence of a flat energy band. At higher excitation power, we 
demonstrate the ability to trigger bosonic condensation 
of polaritons in the flat band. Combined interferometric 
and spectral measurements reveal that condensation oc¬ 
curs on highly localized modes, reflecting the elementary 
eigenstates produced by geometric frustration. These lo¬ 
calized modes arise from disorder in the sample, to which 
flat band states are extremely sensitive due to their infi¬ 
nite effective mass. 

Our ID Lieb lattice of coupled micro-pillars [Fig. Ef] is 
obtained by processing a planar microcavity (of nominal 
Q-factor 70000) grown by molecular beam epitaxy. The 
cavity consists of a A /2 Gao. 05 Alo. 95 As layer surrounded 
by two Gao. 2 Alo. 8 As/Gao. 05 Alo. 95 As Bragg mirrors with 
28 and 40 pairs in the top/bottom mirrors, respectively. 
Twelve GaAs quantum wells of width 7 nm are inserted 
in the structure, resulting in a 15 meV Rabi splitting. 
Micropillars (see close-up in Fig. [l^i) are patterned by 
dry etching down to the GaAs substrate. The diameter of 
each pillar is 3 /im, and the distance between two adjacent 
pillars is 2.4 /im, so that they spatially overlap, allowing 
for the tunneling of polaritons □sin. 

The unit cell of the ID Lieb lattice (Fig. Ef) contains 
three sites (A,B,C) linked by a coupling constant t\ be¬ 
tween A and B, and 1 2 between B and C. In general, this 
lattice exhibits three dispersive bands, as shown in the 
tight-binding calculation of Fig. [1J (dashed lines). In 
the particular case of sites A and C having equal ener¬ 
gies, the middle band is flat and gapped from the two 
remaining dispersive bands (solid lines). The existence 
of the flat band is independent of the energy of pillars B 
and of the couplings £i,£ 2 - 

To probe the polariton properties in this lattice, we ex¬ 
cite it non-resonantly with a continuous-wave monomode 
laser at 740 nm. The spot is elliptical with 2 /im-width 
and 60 /im-length (intensity FWHM), corresponding to 
12 unit cells of the lattice. Experiments are performed at 
6 K and —8 meV cavity-exciton detuning. The emission 
of the sample is collected with a 0.5 numerical aperture 
objective and focused on the entrance slit of a spectrom¬ 
eter coupled to a CCD camera. Imaging of the sample 
surface (resp. the Fourier plane) allows for studying the 
polariton modes in real (resp. reciprocal) space. We re¬ 
solve the emission in polarization, and consider either the 
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FIG. 2. b,h, Energy-resolved far field emission of the struc¬ 
ture at low excitation power, for polarization parallel (TM) 
and perpendicular (TE) to the lattice. Black lines are fits us¬ 
ing the tight-binding Hamiltonian 'Hueb [see Eq. (m)]. When 
the pump excites equally all type of pillars (a), condensation 
occurs in the upper dispersive band and is TM polarized: b-d 
monitor the far field emission when increasing the pump power 
from P = 0.05 P t h to P — 2P t h, where P t h is the threshold 
power (P t h — 4 mW for this condensation process). In d con¬ 
densation is achieved at the top of the upper dispersive band, 
e, Corresponding real space image of the condensate, f, Calcu¬ 
lated trajectory of the eigenvalues of Eq. (|T]) in the complex 
plane, as the power is increased (see arrows), g, When the 
pump is centered on the line of pillars A, condensation occurs 
in the flat band and is TE polarized: h-j monitor the far field 
emission with increasing pump power (P t h — 5 mW for this 
condensation process), k, Corresponding real space image of 
the flat band condensate. 1, Same calculation as f but for the 
pump configuration g, reproducing condensation in the flat 
band. 


polarization parallel (TM) or perpendicular (TE) to the 
lattice (Fig. EL 

Under low pumping power, incoherent relaxation of 
polaritons results in the population of all energy bands, 
allowing a direct visualization of the band structure. 
The corresponding far field emission is shown in Figs. It 
and §1 for emitted light polarized TM and TE, respec- 
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FIG. 3. a-c, Interferograms obtained by superposing two images of the flat band condensate (at P — 2 P t h) shifted by Ax = 0, 
2 and 4 unit cells, d-f, Same experiment realized with the upper dispersive band condensate, g, First-order spatial coherence 
g ^ deduced from the visibility of the interferograms, as a function of Ax. h, Solid lines: Half participation ratio HPR calculated 
from Tihieb as a function of (diagonal) disorder strength, in the flat band and at the top (k = 0) of the upper dispersive band. 
Shaded aeras denote the standard deviation, and the vertical line indicates the estimated disorder strength in our sample. Dots: 
HPR in the flat band calculated from the full Hamiltonian 'H(P) [Eq. Q]. 


tively. In TM polarization a gapped flat band is clearly 
observed, while in TE polarization the middle band is 
dispersive and crosses the upper band. This difference 
arises from polarization-dependent boundary conditions 
for the photonic modes of the pillars, which induce 
a spectral detuning between pillars A and C for TE 
polarization. Both band structures are well reproduced 
(black lines in Figs. [ 2 | 3 ,h) by a tight-binding Hamiltonian 
including next-nearest-neighbor couplings t! \ Hueb = 

E \lj) (l 0 1 - £-(U I A,-) (Bj \+t 2 \Bj) (Cj \+t 2 I Cj) (B j+1 1) 
“EjQlQ) (Bj +ll + !Q) (Q+il + \ A j) (QI + IQ) (Aj+i|)+ h.c. 
where \lj) denotes the wavefunction (of energy E{) of 
pillar l = A, P, C in the j th unit cell. 

The real space pattern corresponding to the flat band 
of Fig. |2j> (i.e. below the condensation threshold) can 
be reconstructed by spectrally filtering the image. The 
result is shown in Fig. EL and compared to the tight- 
binding prediction of Fig. [l^. A characteristic structure 
is observed, with sites B being dark, i.e., containing no 
particles. This is a direct signature of geometric frustra¬ 
tion: a destructive interference between sites A and C 
cancels the net flow of particles to sites B and prevents 
propagation of flat band states [27] [39] . 

We now consider polariton condensation [40] in this ID 
Lieb lattice. Upon increasing the pump power, conden¬ 
sation is triggered when the relaxation towards a given 
mode becomes faster than its decay mt The mode be¬ 
comes macroscopically occupied and spontaneous coher¬ 
ence sets in. This process can be understood by analyzing 
the non-Hermitian operator derived from the generalized 
Gross-Pitaevskii equation [42] : 

U{P) = Uu eb + E [%fijP + | - 7;)] \lj) (h\ 

l=A,B,C 


The second part in Eq. 0 is complex-valued and de¬ 
scribes the effect of pump and dissipation which adds to 
the lattice Hamiltonian Hueb that is fixed by the fabri¬ 


cated structure. The first term accounts for the repulsive 
interaction between polaritons and reservoir excitons, 
with gn being the corresponding interaction constant, 7 r 
and R the exciton decay and relaxation rates, respec¬ 
tively. This repulsive interaction produces a blueshift of 
the pillar energies, proportional to the pump power P; 
fij denotes the fraction of the total pump power on each 
pillar. The second term describes the dissipation rate of 
each pillar: it is the sum of the passive dissipation rate 7 1 
imposed by polariton decay, and a gain term proportional 
fo P. 

The eigenvalues of 'H(P) are complex-valued uj n = 
v n —i^n/2 and describe the frequency and damping of the 
linear fluctuations around the uncondensed state [35]. As 
the power P is ramped up, all complex eigenvalues flow 
towards the real axis from below (see Figs. [2]f,l) . The 
condensation threshold for a given mode is reached when 
its eigenvalue uj n crosses the real axis (net loss y n = 0 ), 
meaning that gain overcomes the decay of polaritons. 
Whereas the polariton decay is fixed by the parameters of 
our structure, the gain is proportional to the spatial over¬ 
lap with the pump and can thus be controlled by tailoring 
the pump profile (described by fij) [35] . This is indeed 
what we demonstrate in Fig. [2] in the ID Lieb lattice 
we can control condensation to occur either in the upper 
dispersive band (left column) or in the flat band (right 
column), by adjusting the pump spatial configuration. 

We first pump the lattice in a symmetric manner, as 
sketched in Fig. [ 2 ^: all types of micropillars (A,B,C) are 
pumped by the same amount (/'Aj—fuj—fcj)^ Figs.^-d 
monitor the momentum space emission in TM polariza¬ 
tion, when increasing the pump power from P = 0.05 Pth 
to P = 2 P t h, where P t h = 4 mW is the threshold power 
(corresponding data for TE polarization are shown in 
the Supplemental Material [43]). Polariton condensation 
eventually occurs in the upper dispersive band (Fig. HU 
with the emission collapsing at the center of the Brillouin 
zones m and showing narrow linewidth. This conden- 
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sation process is well reproduced by calculating the tra¬ 
jectory of the eigenvalues of T-L(P) in the complex plane 
(see Fig. §)• The parameters of 7-Lueb are kept constant 
and the evolution of the eigenvalues with pump power 
is completely driven by the second part of Eq. ([!]). We 
observe that the k = 0 upper band mode (enclosed by 
the square box) indeed reaches the real axis first for this 
pump configuration. 

In the following, we show that asymmetric pumping 
of the same structure allows triggering condensation into 
the flat band. Compared to upper band modes, flat band 
modes have larger amplitude on pillars A (Fig. EM- 
Pumping favorably these pillars ( fs^fcj ) thus 
enhances the gain and lowers the condensation threshold 
of flat band modes, as confirmed by the calculated spec¬ 
trum in Fig. [2j. Using such asymmetric pumping, pillars 
A develop an additional blueshift with respect to pillars 
C, which destroys the flat band in TM polarization [43] . 
Nevertheless, in TE polarization this blueshift compen¬ 
sates the photonic detuning between pillars A and C so 
that a gapped flat band is formed, as experimentally ob¬ 
served in Fig. [2j. Increasing the power further triggers 
condensation in this flat band (Fig. [2]). The emission in 
real space (Fig. [2jx) shows dark B sites characteristic of 
geometric frustration, in clear contrast with the disper¬ 
sive band condensate (Fig. [2^). 

A peculiarity of a flat band is that, due to its macro¬ 
scopic degeneracy, eigenstates of arbitrary lengths can 
be formed. Among these, a set of maximally localized 
eigenmodes can be defined [27], which for the Lieb lat¬ 
tice extend over 3 pillars only [39]. An example of such 
a plaquette-state is highlighted in gray in Fig. [TJd. Note 
that delocalized eigenstates can also be obtained by lin¬ 
ear superposition of different plaquette-states. Thus, the 
question arises of which kind of state is actually picked 
up for condensation in the flat band. 

To answer this question, we use interferometric mea¬ 
surements and high spectral resolution imaging. We first 
investigate the first-order spatial coherence g^(Ax) of 
the flat band emission, by superposing two images shifted 
by Ax along the periodic direction, and extracting the 
averaged visibility of the resulting interference fringes 
m- Figs. §l-C show the measured interferograms for 
increasing values of Ax: the visibility decays exponen¬ 
tially (Fig. [3^) with a characteristic length of 1.6 db 0.3 
unit cells. This indicates the presence of many indepen¬ 
dent condensates, localized on 1.6 unit cells on average, 
i.e. close to plaquett e-states. To confirm this picture we 
use high spectral resolution imaging. Fig. & displays the 
energy-resolved emission of the flat band condensate in 
real space, along the line of pillars A. By repeating the 
measurement for the line of pillars B and C [43] 5 we color- 
code in Fig. [ 4 J 3 the emission energy measured on each 
pillar. It shows significant spectral variations along the 
lattice, allowing to visualize the localized condensates. 
A general measure of the localization of these conden¬ 
sates is provided by the half participation ratio, defined 


as HPR = \’ t Pcj\ 2 ) 2 /(2Y,j IV’Cjl 4 ), where tp C j is the 

wavefunction amplitude on site C of the j th unit cell. The 
HPR quantifies the portion of the lattice where the wave- 
function differs markedly from zero, and coincides with 
the localization length in the case of exponentially local¬ 
ized states 02- I n Fig. starting from the left, one 
successively identifies flat band condensates with an HPR 
of ~ 1.5, 1 and 2 unit cells. This confirms that flat band 
condensation is multimode, occurring simultaneously on 
strongly localized and independent modes. 



FIG. 4. a, Energy-resolved emission of the flat band conden¬ 
sate in real space (see Fig. |2jc) along the line of pillars A. We 
colour-code in b the emission energy along the lattice: several 
localized modes are visible, emitting at distinct energies, c, 
Energy-resolved emission of the upper dispersive band con¬ 
densate in real space (see Fig. [2^), along the line of pillars B 
and C. d, Corresponding spatial map of the emission energy, 
showing monochromatic emission over the whole condensate. 

Let us compare these observations to the situation 
where condensation takes place in the upper dispersive 
band. Keeping the very same location on the sam¬ 
ple, Figs. -f show the corresponding interferograms ob¬ 
tained when the system is pumped symmetrically (Fig. 

2 r). The deduced decays much more slowly (Fig. 

3 5 ), revealing highly extended spatial coherence. Con¬ 
sistently, the energy-resolved measurements of Fig. B-d 
show that the emission is monochromatic over the whole 
condensate within the experimental resolution (30 geV): 
condensation here occurs on a single extended mode, hav¬ 
ing an HPR > 4 unit cells. 

This striking difference between both types of conden¬ 
sates can be understood by taking system disorder into 
account. A flat band is robust to non-diagonal disorder 
(i.e. on the tunnel couplings) [T9] , but any small amount 
of diagonal disorder (on the on-site energies) breaks the 
macroscopic degeneracy: extended eigenstates cannot be 
formed anymore, resulting in a strong localization effect 
m nn Eg. This can be seen by calculating the eigen¬ 
states of the lattice Hueb in presence of diagonal disor¬ 
der. We plot in Fig. It the corresponding half partici¬ 
pation ratio HPR. In the flat band (red line) the HPR 
is very small and essentially independent of the disorder 
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strength, illustrating the non-perturbative effect of dis¬ 
order. Its value HPR ~ 1.8 =b 1 unit cells reproduces 
well the experimental HPR. By contrast, the HPR in 
the upper dispersive band (black line) decays as a power 
law with increasing disorder, and is several times larger 
than in the flat band for weak disorder. We calculate 
HPR = 8.5 =b 4 for the estimated disorder strength of 
our sample (30 /xeV, vertical line in Fig. [3jn) , which is 
compatible with the extended condensate experimentally 
observed. 

To investigate the possible influence of pump and dis¬ 
sipation on the strong localization observed in the flat 
band, we also calculate the HPR using the complete 
Hamiltonian H(P) [Eq. 0], In absence of disorder, we 
find that the HPR in the flat band is bounded only by 
the finite size of the pump spot [45] (as is the case for the 
upper band), showing that pump and dissipation alone 
cannot lead to the observed localization. When includ¬ 
ing disorder, the HPR (red dots in Fig. 0) is close to 
the calculation based on %Lieb only. This allows con¬ 
cluding that pump and dissipation play a minor role in 
the observed localization effect. We note that polariton- 
polariton interactions (cubic nonlinearity) could also lead 
to localization in the flat band [23lf25j , but in the current 
experiment this nonlinearity is of order ^ 0.1 fieV and 
thus much smaller than disorder (30 /xeV) [43] . We thus 
neglected this term in the analysis. 

Finally, the experimental reproducibility of our results 
has been verified on ten different lattice realizations. 
Each time the disorder landscape, which allows the for¬ 
mation of an extended monomode condensate in the up¬ 
per band, leads to a fragmentation of the flat band con¬ 
densate into highly localized modes. To conclude, our 
experiments demonstrate the extreme sensitivity of a flat 
band condensate to diagonal disorder, due to the infi¬ 
nite effective mass in such band. Bosonic condensation 
in a flat energy band opens up novel perspectives in the 
simulation of mesoscopic and condensed matter phenom¬ 
ena. In particular, access to many-body flat band physics 
123 EE Hz] could be gained using resonant pumping, 
where the interaction energy is set by the spectral de¬ 
tuning between the pump and the considered quantum 
state [48]. Furthermore, spin-orbit coupling has recently 
been engineered in polaritonic systems [49]. This should 
allow to further explore the delicate interplay of frustra¬ 
tion, interactions and topology [50j EB in macroscopic 
quantum states. 
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